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Abstract 

We derive three-dimensional integrable mappings which have two invariants. 

1 Introduction 

In this paper we focus on three-dimensional integrable autonomous mappings preserving at least one 
three-quadratic (possibly rational) integral (we considered the four-dimensional case (together with 
its generalization to higher dimensions) in [S|). A major reason for such a study is the lack off 
results on three-dimensional integrable mappings. A recent paper, which makes some progress on 
three-dimensional integrable mappings, is |Sj • In this paper Hirota et al used algebraic entropy 0] to 
determine which three-dimensional mappings of a particular form had polynomial growth implying 
zero algebraic entropy. Having discovered all such possible mappings, they used a procedure outlined 
in their paper to find two functionally independent conserved quantities for each map. In this paper 
we will take a different approach to the one used by Hirota et al to construct three-dimensional 
integrable mappings. For the purposes of this paper, we consider a three-dimensional autonomous 
mapping integrable if there exist two functionally independent integrals in involution with respect to 
some Poisson structure. 

The plan of this paper is as follows : In Section 2 we derive three-dimensional volume-preserving 
mappings which preserve a three-quadratic expression (a method introduced in [2] on a rational four- 
quadratic expression), then assuming that these three-dimensional volume- preserving mappings have 
a second integral with a particular ansatz we find 3 three-dimensional volume-preserving integrable 
mappings. In Section 3 we use the processes of reparametrization and replacement |S1 |5j (terms 
introduced and defined in [5]) to construct three-dimensional measure-preserving integrable mappings. 

2 Three-Dimensional Volume-Preserving Mappings 

In this section we construct three-dimensional volume-preserving mappings (orientation-reversing and 
-preserving) 1 possessing two integrals, at least one of the integrals being quadratic in the three vari- 
ables. 

We begin with the orientation-reversing case. Consider the three-quadratic expression 

I{x,y,z) = Y J A ai , 3lll Z 1 / 1 ^ 1 , (01,01,71 = 0,1,2), (1) 

where A ai/ 3 l7l are independent parameters. Assume that Q is invariant under a cyclic permutation 
of variables 2 , i.e. I(x,y,z) — I(y,z,x), and that the mapping, L, preserving I(x, y, z) is reversible, 
i.e. 

L o G o L = G, (2) 
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1 A mappings L is orientation reversing (orientation preserving) if det dL = — 1 (detdL = 1). 

2 This guarantees that the mapping preserving this integral takes the form x' = y, y' = z, z 1 = F(x, y, z), where F is 
some function. 
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with reversing symmetry G : x' = y,y' — x,z' = x. The reversing symmetry also implies that 
I{x, y, z) — I(z, y, x). Under these conditions we obtain the integral 

I(x,y,z) = Aix 2 y 2 z 2 + A 2 xyz(xy + xy + yz) + A^(x 2 y 2 + x 2 z 2 + y 2 z 2 ) 

+ A±xyz(x + y + z) + A§(x 2 y + x 2 z + xy 2 + xz 2 + y 2 z + yz 2 ) 

+ A 6 xyz + A 7 (x 2 + y 2 + z 2 ) + A 8 (xy + xz + yz) + A 9 (x + y + z). (3) 

The mapping, L, which leaves the integral © invariant can be derived by setting x' — x and y' = y 
and differencing the integral ©, i.e. I(x' ,y' , z') = I(x,y,z') = I(x,y,z). Then assuming z' ^ z, we 
can solve for z 1 , to obtain the involution L z . Finally, composing L z with the cyclic shift L c : x' = 
y, y' = z, z' = x, i.e. L = L z o L Cl we obtain the non-trivial volume-preserving orientation-reversing 
mapping, L, 

x' = y 

y = z 

, _ A 2 y 2 z 2 + A 4 yz(y + z) + A 5 (y 2 + z 2 ) + A 6 yz + A 8 {y + z) + A 9 



A iy 2 z 2 + A 2 yz(y + z) + A 3 {y 2 + z 2 ) + A 4 yz + A 5 {y + z) + A 7 ' 

For a slightly more general map see 

Next, assume that the mapping (0} has a second integral with the following ansatz 3 



(4) 



I(x,y,z) = J2 A «^ x^y^z 72 , (a 2 , 72 -0,1,2, f3 2 = 0,1,2,3,4). (5) 

where A a2 p 2l2 are independent parameters. As the mapping (0J is reversible we also have I%(x, y, z) = 
I 2 (z, y, x). We have found the following mappings which simultaneously preserve integrals of the form 
© and © : 

x' = y 

y' = z (6) 

z , = _ x _ P{y + z) 2 + e{y + z) + j 
p(y + z)+ 1 

with integrals 

I\ = P{x 2 y + x 2 z + xy 2 + xz 2 + y 2 z + yz 2 + 2xyz) + ^(x 2 + y 2 + z 2 ) 

+ e{xy + xz + yz) + £(x + y + z) (7) 

h = [3 2 (x + y) 2 {y + z) 2 +[3(e- 1 )[(x + y) 2 (y + z) + {x + y)(y + z) 2 } 

+ 7 (e - 2 7 ) [(a: + yf + {y + zf\ + + (e - 7 )(e - 2j)](x + y)(y + z) 

+ Z(e-2 1 )[(x + y) + (y + z)} (8) 



and 



x = y 

y> = z (9) 
2f3yz + e(y + z) + £ 



z = —x 



ayz + (3(y + z) + 7 
with integrals 

11 = a 2 x 2 y 2 z 2 + 2af3xyz{xy + xz + yz) + aexyz(x + y + z) 

+ P 2 (xy + xz + yz) 2 + (ie{x 2 y + x 2 z + xy 2 + xz 2 + y 2 z + yz 2 + Axyz) 

+ (a£ — 2(3j)xyz + 7 e(a: 2 + y 2 + z 2 — xy — xz — yz) — "f 2 {x 2 + y 2 + z 2 ) 

+ {fc + £ 2 )(xy + xz + yz) + C(e - j)(x + y + z) (10) 

1 2 = {ft + e)[axyz(x — y + z) + j(x 2 + y 2 + z 2 — xy — yz + xz) + £(x + z)] 

+ [I3 2 (x + z)+ f3e(x + z + 1) + ( 2 ][xy + xz + yz - y 2 }. (11) 
3 See the Appendix for the reason why the ansatz has this form. 
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We next consider the orientation-preserving case. Consider the three-quadratic expression £[J 
possessing the symmetry I{x, y, z) — I(y, z, —x) 4 . Following the procedure outlined above we obtain 
the mapping 



y' = z 



V (12) 

x A{y-z) 
+ Byz + C 

with integrals 

l x = B 2 x 2 y 2 z 2 - ABxyz(x - y + z) - A(A + C) (xy - xz + yz) 

-C(A + C)(x 2 +y 2 + z 2 ) (13) 
h = B a x 2 y 2 z 2 -C 3 (xy + xz + yz) + ABC{2xz-2xy 2 z-y) + A 2 (B-C)y 2 

- BC 2 (x 2 + y 2 + z 2 + x + y + z + xyz[x + y + z]) 

+ 2AB 2 xy 2 z - B 2 Cxyz(xyz + 1) - AC 2 {x + y){y + z) (14) 
We close this section with the following remark. 

Remark The mapping © under the coordinate transformation X = x + y,Y = y + z can be 
reduced to a two-dimensional area-preserving mapping, i.e. 

L X : X' = Y, Y> = -X- ie -^ + Z (15) 

with the second integral, I2, becoming 

I = /3 2 X 2 Y 2 + P(e--/)(X 2 Y + XY 2 )+ 1 (t-2 1 ){X 2 + Y 2 ) 

+ m + (e - 7)(e " 2l)]XY + £(e - 2j)(X + Y). (16) 

Note that the first integral, 7i, does not reduce under this transformation. 

In fact, the mapping © is a member of a recently-discovered hierarchy of integrable mappings 
given in [H], the three-dimensional asymmetric mapping 5 being 

/% + z) 2 + e(y + z) + g 

x = — x — 

(3{y + z) + 70 

p(x' + zf + e{x' + z) + Ci 



&{x' + z) + 71 
, _ _ f3(x' + y') 2 + e(x' + i/)+^ 

P{x>+y') + l2 ■ (U) 

3 Three-Dimensional Measure-Preserving Mappings 

In this section we apply the processes of reparametrization and replacement to the three-dimensional 
volume-preserving integrable mappings constructed above to construct measure-preserving integrable 
mappings. These examples illustrate how integrable three-dimensional families of mappings can be 
embedded in larger (i.e. higher number of parameters) integrable three-dimensional families of map- 
pings via the processes of reparametrization and replacement. 
Consider the mapping © when a = and e = 7, i.e. 

x' = y 
/% + z) + 7 



4 This symmetry is due to 

5 In |B| this asymmetric mapping was shown to be obtained as a composition of three involutions. We believe that 
this guarantees the reversibility of this mapping. It seems, more generally, that a mapping obtained from a composition 
of involutions is reversible, see [B] for examples of such mappings. 
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which has integrals 

I\ = 0(xy + xz + yz) 2 + 7(2; + y)(x + z)(y + z) + £(xy + xz + yz) (19) 
h = (0 + j)[-f(x 2 + y 2 + z 2 — xy - yz + xz) + £(x + z)] 

+ [I3 2 {x + z)+ 0j(x + z + 1) + -f 2 ][xy + xz + yz - y 2 }. (20) 

Notice that the parameters 0, 7 and £ now appear linearly in the integral I±. Reparametrizing the 
parameters, i.e. — » 0o+0iK,~/ — > 7o + 7i-^>£ £,q+£,iK and the integral I\ — > Ji = Ji +/io + Mi^ 5 
we obtain the mapping 



y 



y' = 



, = 2(/3 + (3iK)yz + (70 + liK)(y + z) + + gjjf 

(/3 + /3iif)(y + z) + (7o+7i^) ' 

Using I\(x,y,z) = (as I\(x,y,z) = =>• Ix(x',y', z') = 0) a new integral if = k(x,y,z) can be 
defined. Define the map Ljf to be the map (|21[1 with replacement K — k(x,y,z). The map Lk has 
two integrals k(x,y,z) and 7 2 = h(x,y,z)\K=k{x,y,z)> i- e - 

/3o(a:T/ + + yz) 2 + 7o(x + y)(ai + z)(y + z) + £ (xy + xz + yz) + fi Q 
k = — - 



0i (xy + xz + yz) 2 + 71 (x + y) (x + z) (y + z) + £1 (xy + xz + yz) + fii 



(22) 



h = {[Po + 7o + (f3i+li)K}[{"f + JiK)(x 2 +y 2 + z 2 -xy-yz + xz) 
+ (£0 + ZiK)(x + z)] + [(0o + 0iK) 2 (x + z) 
+ (0 O + 0iK)(r/o + nK){x + z + 1) 

+ (7o + liK) 2 ][xy + xz + yz- y 2 ]}\ K =k(x, v ,z)- (23) 
The map Lk is also measure preserving with 

m(x,y,z) 



(24) 



dK 

Consider the mapping when a — and e = 7 = 0, i.e. 

x' = y 

z> = _ x _mvz + v + z) + t (25) 
0{y + z + i) K ' 

which has integrals 

Ix — 0[(xy + xz + yz) 2 + (x + y)(x + z)(y + z)] + £(xy + xz + yz) (26) 
h = P(x + z)(x + z + xy + xz + yz-y 2 )+£(x + z). (27) 

Notice that the parameters and £ now appear linearly in both integrals. Reparametrizing the 
parameters and the integrals, i.e. — > 0q + 0\K\ + P2K2, £ —> Co + + £2-^21 h ~* h — 

h + Mo + ViK 1 + ^2^2 and I 2 — > h = h + ^0 + V1K1 + v 2 K 2 , we obtain the mapping 

x' = y 

y' = z 

, = (00 + 0xKx + 02K 2 )(2yz + y + z) + Co + giifi + &K 2 

(0 o +0lK 1 +0 2 K 2 )(y + z + l) ' 1 8j 

with integrals h(x, y, z) and hix, y, z). Setting ii(a;, y, z) = and ^(a:, y, z) = it is possible to solve 
for Ki = ki(x, y, z) and K 2 = k 2 (x, y, z) as I\ and J 2 are linear in ^ and K 2 . Define the map Lx 1 k 2 
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to be the map (|2"5)l with replacements K\ = ki(x, y, z) and K 2 = k 2 (x, y, z). The map Lk x k 2 has the 
integrals Ki — ki(x,y,z) and K 2 — k 2 {x,y, z). The map Lx x k 2 is also measure preserving with 



m(x,y,z) 



dh 
dKi 
dh 



dh 
dh 

OK; 



(29) 



The integrals to the above maps can be shown to be functionally independent and in involution with 
respect to the following Poisson structure 



/ |£ -|£ 



m(x,y, 



\ 



_ dl_ 

dz 
01 




(30) 



where / is either one of integrals and m(x, y, z) is the measure. 

Finally, we consider the mapping JBJ . As noted in the remark at the end of Section 2 we can use a 
coordinate transformation, i.e. X = x + y and Y = y + z, to reduce the mapping to a two-dimensional 
mapping. Importantly, however, the three-quadratic integral, I\, does not reduce under this coordinate 
transformation and as a result if we use the processes of reparametrisation and replacement on this 
integral then the resulting mapping, L r , is not reducible to a two dimensional mapping, although for 
every fixed K it is. The remark at the end of Section 2 also shows that the reduced mapping has a 
biquadratic integral and thus can be explicitly integrated 6 , see [Jj. This result can be used to integrate 
the mapping L r also but this time curve- wise (leaf- wise). 7 
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Appendix 

Our work on multidimensional integrable mappings (particularly 6 ) has lead us to the following ob- 
servation about the degree of the variables that occur in the integrals : 

A 2n- dimensional volume-preserving integrable mapping with variables vq, . . . , v 2n -\, which has 
one n-quadratic integral, possesses an additional (n — 1) integrals of the following form 



h = J2 A a 21 ...a 2:2n VQ 21 . . .«2n-l > ( a 21, ■ ■ . , £*2,2n = 0, . . . , 4) 



In = Y!,Aa n i...an,2n V n ' 1 • • • V 2n- 1 i ( a n,l, ■ ■ ■ ,Qn,2n = 0, . . . ,2n). 



(31) 



While a (2n+l)- dimensional volume-preserving integrable mapping with variables Vo, . . . ,v 2n , which 
has one n-quadratic integral, possesses an additional n integrals of the following form 

h = E^W...a 2 , 2 „ + 1 Vo 21 ' ■ - V 2n 2n+1 . ("21, ■ • ■ , Oi 2 . 2n+1 = 0, . . . , 4) 

: (32) 

In+l = J2 ^a„ + i,i...Q„ +li2 „ + i"o" +1,1 ■ ■ ■ V 2 " +1 ' 2 ™ + 1 , (ttn+1,1, ■ ■ • , a n +l,2n+l = 0, . . . , 2n + 2) . 

In the case we have considered in this paper, the power of the first and last variables ranges from 
to 2. 



6 This also is true for its asymmetric form. 

7 We believe that the maps considered in Case 2 of |S] Section 3] can be integrated in this way also. 
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